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ABSTRACT 

Binary neutron stars merger (NS^M) at cosmological distances is probably 
the only 7-ray bursts model based on an independently observed phenomenon 
which is known to be taking place at a comparable rate. We describe this model, 
its predictions and some open questions. 

Cosmological 7-Ray Bursts and Fireballs 



Compton-GRO has demonstrated, quite convincingly, that 7-ray bursts 
(grbs) originate from cosmological sources^ '^. Evidence for the predicted'^'^ cor- 
relations between the duration, the strength and the hardness of the bursts 
begins to emerge^'^. Preliminary analysis suggests that the weakest bursts orig- 
inate from z ?a 1, in agreement with a cosmological C jCrain distribution^, corre- 
O ! sponding to a local rate of ~ 10~^/year/galaxy (depending on the cosmological 

CQ . model and on other factors) . The energy released in each burst depends also on 

the cosmological model 10^'^>-E<10^^ergs if the energy emission is isotropic. 
Q ■ The intense energy released in a small volume (evident by the rapid rise 

Tj- \ time of some of the pulses) implies that any cosmological grb source is initially 

^ \ optically thick^ to 77 —^ e^e~ . The large initial optical depth prevent us from 

^ . observing directly the photons released by the source regardless of the specific 

Qh| nature or the source. The sources produce an optically thick radiation-electron- 

Q . positrons plasma "fireball", which behaves like a fluid, expands and reaches 

^ \ relativistic velocities''^. The observed radiation emerges only after the fireball 

^ ■ has expanded significantly and became optically thin. 

One should divide, therefore, the discussion of cosmological grbs to a 
discussion of the nature of the energy source (for which we present a model 
'k>( \ here) and a discussion of the fireball phase (which we address elsewhere in this 

'^ ■ volume). For the paper it is sufficient to recall that the fireball must reach 

^ ' ultra-relativistic velocities with a Lorentz factor 7>10^ to produce a grb. Since 

7 Ri E IM(? (where E is the total energy of the fireball and M is the mass of 
the baryons in the fireball) the condition 7 > 10^ sets a strong upper limit on 
the amount of baryons: M < Ej-^c^ ^ .blO'^ Mq{E /10^^ergs){-f /lO'^)-'^ . This 
condition poses a strong constraint on grb models. 

NS^M and GRBs - Agreement at a Glance 

Neutron star binaries, such as the one observed in the famous binary pul- 
sar PSR 1916-M3, end their life in a catastrophic merge event (denoted here 
NS^M). Using the three observed binary pulsars we can estimate the expected 
rate of NS^M events^'^'^^ as ~ 10~^'^^'^/year/galaxy. An energy comparable to 
a neutron star binding energy (>5 x 10^^ ergs) is released in NS^Ms mostly as 
neutrinos and gravitational radiation. The neutrino signal is comparable in its 



signature to supernova neutrino signals which are thousand times more frequent. 
It is unhkely that it wiU ever be detected. The gravitational radiation pulses, 
have however, a unique signature and Two gravitational radiation detectors, 
LIGO and VIRGO are currently constructed to detect them. 

Several years ago Eichler, Livio Piran and Schramm^"^, (see also I4,i5,i6,i7^ 
suggested that grbs originate at NS^Ms. Between 10~^ to 10~^ of the total 
energy released in NS^Ms is sufficient to power a grb at a cosmological dis- 
tance. The required energy could be converted to electromagnetic energy either 
vial3, 15 vv -^ e~^e~ or via magnetic processes in an accretion disk that forms 
in the merger^^. The rates of NS^Ms estimated from binary pulsars and the 
observed rate of grbs measured by BATSE and estimated from cosmological fits 
are within half an order of magnitude from each other. A remarkable agreements 
in view of the large uncertainties involved in both estimate. 

Numerical Simulations of NS^M - Some Answers to Further Questions 

It worthwhile, therefore, to explore whether the mergers can produce clean 
enough fireballs (i.e. fireballs with sufficiently low baryonic load) as required 
from the fireball analysis and to ask whether enough energy can be converted to 
electromagnetic energy in this events. To address these iissues we^^ developed 
a numerical code that follows neutron star binary mergers and calculates the 
thermodynamic conditions of the coalesced binary. The process of coalescence, 
from initial contact to the formation of an axially symmetric object, takes only 
a few orbital periods. Some of the material from the two neutron stars is shed, 
forming a thick disk around the central, coalesced object. The mass of this 
disk depends on the initial neutron star spins; higher spin rates resulting in 
greater mass loss, and thus more massive disks. For spin rates that are most 
likely to be applicable to real systems, the central coalesced object has a mass 
of 2.4Mq, which is tantalizingly close to the maximum mass allowed by any 
neutron star equation of state for an object that is supported in part by rotation. 
Using a realistic nuclear equation of state we estimate the temperatures after 
the coalescence: the central object is at a temperature of ~ lOMeV, whilst the 
disk is heated by shocks to a temperature of 2-4MeV. 

A typical density cut perpendicular to the equatorial plan is shown in 
Fig. 1. The disk is thick, almost toroidal; the material having expanded on 
heating through shocks. This disk surrounds a central object that is somewhat 
fiattened due to its rapid rotation. An almost empty centrifugal funnel forms 
around the rotating axis and there is practically no material above the polar 
caps. This funnel provides a region in which a baryon free radiation-electron- 
position plasma could form^'^. Neutrinos and antineutrinos from the disk and 
form the polar caps would collide and annihilate preferentially in the funnel (the 
energy in the cm. frame is larger when the colliding u and u approach at obtuse 
angle, a condition that easily holds in the funnel). The numerical computations 
do not show any baryons in the funnels. The resolution of our computation is 
insufficient, howeer, to show that the baryonic load in the funnel is as low as 
needed. The neutrinos radiation pressure on polar cap baryons can generate a 
baryonic wind that will load the fiow. Estimates of this effect^^'^^ show that 
it is negligible if the temperature on the polar caps is sufficiently low. The 
estimated temperature from our computations is ~ 2MeV, which is marginal. 
Our temperature estimate is, however, least certain in low temperature regions 
like this. 



If the core does not collapse directly to a black hole it will emit its thermal 
energy as neutrinos. The neutrino flux is sufliciently large that ~ 10~^ to 10~^ 
of it could be converted to electron-positron pairs via vv -^ e+e~ and produce 
a grb. The time scale for the neutrino burst is short enough to accommodate 
even the shortest rise times observed. An additional energy source that could 
power a grb is the accretion of the disk surrounding the central object. This 
energy source can operate on a longer time scale and it takes place regardless 
of the question of whether the central object collapse directly to a black hole or 
not. 

Open Questions and Predictions 

The numerical calculations support earlier suggestions^^ that the energy 
release in anisotropic and that an empty funnel forms around the rotating axis 
of the binary system. The flreball is highly non spherical and it expands along 
the polar axis and forms a jet. This poses an immediate constraint on the model. 
If the width of the jet is Q than we observe grbs only from a fraction 2^~^ of 
NS^Ms. The rates of grbs and NS^Ms agree only if 6*^0.2 (unless the rate of 
NS^Ms is much higher than the current estimates). A condition which at flrst 
glance is satisfled by the funnel seen in Fig. 1. 

The duration and spectra of grbs vary greatly from one burst to another. 
Both are determined by the fireball phase but the source might contribute in 
producing fireballs with different Lorentz factors and different initial durations. 
Within the funnel the baryonic load will vary as a function of the angular position 
leading to varying final Lorentz factors which, in turn, produce bursts with dif- 
ferent durations and spectra. Another source of variability could arise from the 
interplay between the two energy sources in NS^Ms: Neutrino annihilation and 
accretion energy of the disk. These mechanisms would operate on different time 
scale and produce different looking bursts. An additional source of diversity^^ is 
the distinction between systems that collapse directly to a black hole and those 
that undergo a longer rotating core phase. Finally, black hole-neutron star bi- 
naries are predicted to be as common as neutron star binaries^^. A black hole 
neutron star merger^'^ would produces grbs with different characteristics than 
NS^M. 

NS^M events can take place in a variety of host systems including dwarf 
galaxies, or even in the intergalactic space if the neutron star binary is ejected 
from the host galaxy when it forms^^. Hence, unlike other cosmological models 
it is not essential that an optical counter part will be observed in the location 
of grbs^^. A unique prediction of the NS^M model is that grbs should be ac- 
companied by gravitational radiation signals from the final stages of the merger 
and vice versa (the latter is true only up to the anisotropic emission factor dis- 
cussed earlier). This coincidence could prove or disprove this model. It could 
also serve to increase the sensitivity of the gravitational radiation detectors^^. 
Hopefully, this coincidence will be detected and the model will be confirmed 
when gravitational radiation detectors will become operational at the turn of 
the century. 
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Figure Captions 

Fig. 1 Logarithmic density contour lines at the end of the computation of the 
merger. The contours are logarithmic, at intervals of 0.25 dex (from ^^). 



